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ABSTRACT 



Aims. By neglecting sideways expansion of gamma-ray burst (GRB) jets and assuming their half-opening angle distribution, we estimate the 
detectability of orphan optical afterglows. 

Methods. This estimation is carried out by calculating the durations of off-axis optical afterglows whose flux density exceeds a certain 
observational limit. 

Results. We show that the former assumption leads to more detectable orphans, while the latter suppresses the detectability strongly compared 
with the model with half-opening angle 8j — 0.1. We also considered the effects of other parameters, and find that the effects of the ejecta 
energy Ej and post-jet-break temporal index — Q2 are important but that the effects of the electron-energy distribution index p, electron energy 
equipartition factor e e , and environment density n are insignificant. If Ej and Q2 are determined by other methods, one can constrain the 
half-opening angle distribution of jets by observing orphan afterglows. Adopting a set of "standard" parameters, the detectable rate of orphan 
afterglows is about 1.3 x 10 _2 deg _2 yr _1 , if the observed limiting magnitude is 20 in R-band. 



r Key words, gamma rays: bursts - ISM: jets and outflows - radiation mechanisms: non-thermal 



1. Introduction 

Orphan afterglows are defined as afterglows whose gamma- 
ray bursts (GRBs) are not detected, possibly because of the 
Doppler effect for an off-axis observer. If the GRB afterglows 
are modelled perfectly, the observed rate of orphans and GRBs 
can be used to constrain th e beaming factor of GRBs, as first 
proposed by Rhoad^ dl997h . However, as many parameters are 
not determined well for different afterglows, it will be difficult 
to constrain the beaming factor tightly from optical orphans. 
Because late radio afterglows behave isotropic emission, the 
survey of radio aftergl ows may be helpful for estimating the 
beaming factor ( Levins on et alJ2002tlGal-Yam et alJ2006l) . al- 
though one should be careful to rule out radio transients from 

other sources. 

O ther than constraining the beaming factor jRhoadsl 

1 19971 iDalal. Griest & Pruelll2002t llbtani & Panaitescull2002ir 
some authors have focused on investigations of the de- 
tectability of orphan optica l afterglows both theoreti cally 
JNakar. Piran & Granot|l2002UTotani & Panaitescull2002l) and 



2004) survey o f optical tran s ients, but none was identified 
as an orphan. iRvkoff et al . (2005) performed a 1.5-year 
survey (2003 September to 2005 March) of untriggered 
GRB afterglows. Although no orphan afterglow has been 
observed yet, they give the upper limit of the observed rate 
for a certain hmiting magnitude. The surveys are s till going 
on dMalacrino & Atteial 120061 iMalacrino et al1l200fj|) . On the 
theoretical side. iNakar. Piran & G ranot (2002) considered the 
following afterglow model: all jets propagating in a uniform 
medium (ISM) have a constant initial half-opening angle Oj 
and a constant jet energy Ej. After a jet break takes place 
because the hydrodyna mics of a sideways-expansion jet enters 
an exponential regime (Rhoads 1 19991 ISari. Piran & HabertJ 
Il999l) . the temporal index of light curves becomes — p (where 
p is the power-law index of shock-accelerated electrons). 
This decli ne is too steep for most of the observed late 
afterglows! Liang, & Zhang 2005). On the other hand, many 



experimentally (Hudec 2004; Becker et al. 2004; Rvkoff et al. 
2005 ; Malac rino & Atteial \200(t lRa~G reiner & Schwarj 



Wei & Lu 


2000; 


Salmonson 2003; 


^^mar&Granot 


2003 


Granot & Kumar 2003; Cannizzo. Gehrels & Vishniad 12004) 



20061). iBecker et alJ J20041) give the results of a 5-year (1999- 



show that the sideways expansion of jets is insignificant at the 
relativistic stage. Thus, we consider relativistic jets without 
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Fig. 1. Sketch of afterglow light curves from jets without side- 
ways expansion for an on-axis observer. Dotted, dot-dashed 
and solid lines correspond to three jet opening angles 0.03, 
0.1 and 0.3 respectively, with the same total kinetic energy 
Ej = 1 x 10 51 ergs. The dashed line is the connection of jet 
breaks. 



sideways expansion, and their afterglow light curves for an 
on-axis observer are shown in Fig. ^ The light curves are 
shown in the spherical case and the flux density /„ oc t~ ai 
when the Lorentz factor 7 > 1/6 j. After the jet break time, 
the light curves steepen as f v oc t~ a2 (ai > oti) because of 
the edge effect! Me zsaros &Rees 1999), which is flatter than 
the sideways-expansion case. This will lead to more detectable 
orphan afterglows. A relationship between the jet break time 
and flu x density {f V j oc t~ p ) was found bv lWu. Dai & Liangl 
(2004) analytically and statistically. 

Fig. Q] is plotted based on the fact that the ki- 
netic energies of all GRB jets have a similar value 
jFrail et alJ l200lt iBloom. Frail & Kulkarnil l2003h . The sta- 
tistical standard ene rgy of jets has been discussed by 
several autho rs, e.g., |L ipunov, Postnov & Prokhorov (2001), 
and Panaitescu & Kumar ( 2002JK who conclude that the 



collimation-corrected gamma-ray en ergy has a relatively nar- 
row di stribution, around 5 x 10 50 ergs. lBerger. Kulkarni & Frail 
(2003) also obtained a standard kinetic energy reservoir of af- 
terglows from the statistics on X-ray luminosity. 

Recently, the distribution of the half-opening angle or 
viewing angle was investigated based on several structured-jet 
models JPerna. Sari & Frap 120031: lliang. Wu & Dai l2004 
Nakar. Cranot A Guettal l2004 lOuetta. Piran A Waxmanl 
2005 |. Considering a uniform, sh arp-edge jet (favored 
by lLamb. Donaghv & Grazianil 120051) . whose light curves 
are similar to those in the uni versal structured jet mod- 
els jRossi. Lazzati & Reesl 120021) : and using the observed 
distributi on of half-opening an gles of the jets given 
by lLamb. Donaghv & Grazianil J2005I) . one can derive 
the int rinsic distribut i on of 6j, i.e., P(6j) oc 6*7 1 (see 
also IXu. Wu & Pall |2005). We use this distribution as 



a weight for different openi ng angles as suggested by 
iGuettaTpiran & Waxmar] J2005I) . 

By considering both the effects of the constant half-opening 
angle during jet propagation and the distribution of initial jet 
half-opening angles, we here estimate the detectability of or- 
phan afterglows and find that our results are different from 
the ones in earlier works (e.g. iNak ar. Piran & Granot 2002; 
iTotani & Panaitescull2002h . We present the theoretical model 
in ||2]and give the results of the detectability in ||3] We sum- 
marize our findings and present a brief discussion in f0] 

2. Theoretical analysis 

We consider an adiabatic jet with a total kinetic energy Ej 
and a half-opening angle 8j and neglect sideways expansion. 
The hydrodynamics of the jet behaves as a spherical case 
JSari. Piran & Naravanl [l998). The Lorentz factor of the jet is 
given by 



7(<e) = 8.9(l + ^« D - 1/ '^ 



-l/4 + -3/8 



(l) 



where z is the redshift of the GRB, n the number density of the 
interstellar medium (ISM), and t® d the observed time in units 
of days. We adopt the conventional notation Q = Qk x 10 fc in 
this paper except for special explanations. 

For an on-axis observer, ther e is a break in the l ight curve 
because of the edge effect (Me zsaros & Reesl fT999) when the 
bulk Lorentz factor 7 equals 6 J . The jet-break time is given 
by 

tj = 0.82(1 + z)E) / ^ 1 n~ 1,z el_ 1 days. (2) 

At the jet-break ti me, the flux density in th e slow-cooling case 
{y m < v < v c ) is JWu. Dai & Lian J2004 

F»* = x 50.2 2 - 2 ^ K/ 4r 1)/2 <-i 



V >+1)/Vp-1 (3-p)/12 p(p+3)/3 n -2 
Xe B,-3 S1/6 n ' < 1 



L,28 



x(l+z)(f+ 3 )/ 2 ( — )-(P-D/2 ^Jy, 



(3) 



where « m = 0.73(p - 0.67), K f = 0.09(p + 0.14), and 
k c = (y — 0.46) exp (3. 16 — 1.16p) are the correction fac- 
tors (Granot & Sari 2002J); e e and are the energy equipar- 
tition factors of the electrons and magnetic field, respectively; 
C1/6 — 6(p — 2)/(p — 1); Dl is the luminosity distance; and 
vr = 4.55 x 10 14 Hz is the i?-band frequency taken as the 
observed frequency. On the other han d, in the fast co oling case 
{y c < v), the flux density is dWu. Dai & Lian g 2004) 



F„j = 3508*7^ x SO^-PKfK^^D 



2 

L,28 



p-1 Jp-ij/i^p-l p O+2)/3 _(p+2)/12 
/6 ^j,51 n 



X£ e-l e B,-3 Si/6 ^ 



x(i + z)(f+ 2 )/ 2 (i + r,)- 1 (— )- 



■p/2 



(4) 



where Yj = Y(tj) = + 477e P /e B )/2 is the Com pton 
parameter (Panait escu & Kumarll200o|: Sari & Es"inll200ll) . and 
Tj is the radiation efficiency of electrons. Even though most 
GRB afterglows m atch th e slow cooling case in the statistics by 
IWu~ Dai & Liang (2004), we here consider both fast and slow 
cooling cases, which are different from lNakar. Piran & Granoil 



Y. C. Zou, X. F. Wu & Z. G. Dai: detectability of orphan afterglows 



3 



who only considered the case v > v c > v m for sim- 
plicity. 

For an on-axis observer, the Lorentz factor of the jet 7(t® ) 
at earlier times (t§ < tj) is greater than 9J , so the emission 
properties are the same as those from an isotropic fireball. The 
temporal decay index of the flux density (t@ ) is (2 — 3p) /4 
in th e fast cooling case and 3(1 — p) /4 in the slow cooling 
one rtSari. Piran & Naravanll 998). As index p is mainly in the 
range of 2.0 ~ 2.4, we find that the range of the temporal 
index is about —0.7 to —1.3 for both cases, which is set to be 
a parameter — oc\. When t® > tj and if 9j <C 1 and 7 ;» 1, 
the on-axis observer can only detect a fraction OjJ 2 of the flux 

density in the isotropic fireball case. As j(t@) oc £ ffi 3 ^ 8 [see 
Eq. |0], the late decay index of the flux density a 2 = ai+3/4. 

For an off-axis observer with observing angle 9 obs , the time 
and frequency from on-axis (to, vq) and off-axis (t, v) jets sat- 
isfy to/t ~ v/v q = (1 — /3)/(l — /3cos0 o bs) = a, where 
j3 — y/1 — I/7 2 is the velocity in units of c; thus the flux den- 
sity is 

F u (9 obs ,t) = a?F v/a (0,at), (5) 

in the poin t source approxim ation, which is good enough when 
6 obs > 9j JGranot et al.l2002h . 

Given a limiting flux density f Vl um (corresponding to a 
limiting magnitude mum) for an instrument with a fixed 
exposure time, we can calculate the detectable duration of 
an orphan afterglow: t obs (z,6 obs) 9j,mi im ) = t max - t min , 
where t m i n and t max represent the earlier and later times 
when F v fi obs = F Vt n m . If the maximum observed flux den- 
sity F Vj e oU ,max < Fv,nm, we take t obs = 0. Fig. shows the 
light curves of afterg lows for different observ ing angles. 

Following iNakar. Piran & Granotl ll2002l) . we assume that 
the GRB rate n(z) is proportional to the star formation 
rate (SFR), but we use a different SFR model as follows 
JPorciani & Madaul200ll) . 



n(z) = B- 



JiAz 



Vr> m (i + 2) 3 + r> A 



3.4z 



(6) 



, 22 (1 + z) 1 - 5 

where Q m and Vl\ are the cosmological parameters, and B 
the normalization factor. The parameter B satisfies Rqrb = 
fbR°GRB = Io°(dV/dz)n(z)/(l + z)dz, where f b = 

J2;r t 1 - cos s j)~ 1p vj)Mj/ jSt p ( d ^ is the mean 

beaming factor of GRBs and Rq^ b — 667yr _1 is the observed 
GRB rate. Here we assume all the GRBs can be observed if the 
observer is located within the solid angles of the jets and the 
redshift range is < z < 10. 

If the exposure time is not too long (shorter than t obs ), the 
number of detectable orphan afterglows in a single snapshot 
over the whole sky can be expressed as 
-10 



N, 



orph 



n(z) dV(z) 
(1 + z) dz 



P(9 j )d9 j 



t ob s(z, obs , &j, mii m )9 h s d9 obs (7) 




Fig. 2. Sketch for observations of orphan R-band afterglows. 
The three solid lines indicate the light curves with different ob- 
serving angles. The one with 9 obs = is not an orphan af- 
terglow, which is plotted here as a reference. The horizontal 
dashed line represents a given limiting magnitude. The earliest 
and latest times (i m ; n and i max ) at which an orphan afterglow 
are observed are represented by vertical dotted lines. 



with the upper and lower Umits 9j 



fies ( lLamb. Donaghv & Grazi ani 2005) 



and 6j m i n , which satis- 



ln(^j,max/ ^j,min) 



(8) 



3. Numerical results 

If the model parameters are given, the detectability of or- 
phan afterglows can be estimated by Eq. Q. The main 
difference in detectability comes from the limiting magni- 
tudes of detectors. Fig. [5] shows the number of orphan af- 
terglows that can be detected by one exposure on the whole 
sky. The solid line is the standard result, with parameters 
Ej = 1 X 10 51 erg,n = 1cm- 3 , p = 2.2, a 2 = 1.8, e e = 
0.1, e B = 0.01,1/ = 4.55 x 10 14 Hz, 9 jMn = 0.01, and 
9j may = 1, where Ej , n, y, e P , e r , and v are the same as in 
Nakar. Piran & Granotl J2002h . As the pre -break temporal in- 
dex o f the optical Ught curve is about — 1 ( Zhang & Mezsaros 
|2004 . we choo se a 2 = «i + 3/4 ~ 1.8. The 6» min = 0.01 



max (z, mii m ) is the maximum observing angle, which 



where 9 
satisfies i max 

vational distribution function of half-opening angles of the jets flux density after the break time tj is F v ) ~ F u j (t© / tj ) ~ p 







), and P{9j) is the obser- 



is adopted from Lamb, Donaghv & Graziani (2005). We take 
#max = 1) which does not influence the estimation signifi- 
cantly when we consider the distributio n of half-opening an gles 
of the jets. We also show the results of Nakar. Piran & Granotl 
d2002h in this figure with their canonical ( 9j — 0.1) and op- 
timistic ( 9j = 0.05) parameters. For comparison, we plot the 
dotted line for a fixed half-opening angle 9j — 0.1 and the 
SFR model in INakar. Piran & Granotl J2002I) . One difference 
between the dotted line and the thick dashed line is the temporal 
index after the break time. We can see that the difference in de- 
tectability is due to the sideways expansion. Approximately, the 



4 



Y. C. Zou, X. F. Wu & Z. G. Dai: detectability of orphan afterglows 



Table 1. The ratio of the number of orphan afterglows to the 
total number of afterglows for different limiting magnitudes of 
detectors. The values without and with brackets correspond to 
models A and B in Na kar. Piran & Granot (2002) respectively 
to compare with their results. 



M (MLiuniUkk'i 



■^pcak 


0j 


mi im =23 


mu m =25 


mum=27 


l 


0.05 


0.56(0.68) 


0.73(0.78) 


0.84(0.86) 


l 


0.10 


0.33(0.39) 


0.40(0.58) 


0.63(0.71) 


l 


0.15 


0.26(0.29) 


0.34(0.40) 


0.41(0.59) 


2 


0.10 


0.22(0.36) 


0.34(0.55) 


0.60(0.70) 



in the sideways expansion case and F v (t§) ~ F^,j(t^/tj)~ a2 
in the non-sideways expansion case. Note that F y _j and tj 
have the same values in both cases, since the breaks both 
take place when 7 ~ 1/0 j, and before the j et break time- 
both cases show isotrop ic evolutional behavior (Rhoads 1999; 
iMezsaros & Ree Jl999l) . Neglecting t nlin and 6j in Eq. (7Jl and 
the potential influence of different spectra, we obtain the ra- 
tio of the detectabilities in the two cases (i.e., no sideways 
expansion vs. sideways expansion): N orp h,NSE/N orp h,sE — 
(^,j/^,iim) 11/(8a2) " 3/(2p) - In general, if a 2 < p, then 
ll/(8a2) > 3/(2p) and thus N orp , hNSE > N orph , SE . For 
a larger limiting magnitude (i.e. smaller F Vt n m ), the ratio be- 
comes higher. This is why in Fig. [5] the dotted line is higher 
than the thick dashed line for greater mu m . To show the ef- 
fect of different SFR models, the dot-dashed line uses the SFR 
model in Nak ar. Piran & Granotl J2002I) 1 with z pca k = 1, and 
the dotted line considers the SFR model in Eq. These two 
lines are close to each other, which shows that the effect of the 
SFR models is insignificant. We note that the distribution of 
half-opening angles of the jets leads to further suppression of 
the detectability, which results in the difference between the 
dot-dashed and the solid lines. Combining the effects of the 
sideways expansion and distribution of the jet's half-opening 

angles, we obtain standard results (soli d line in Fig.[3l>. 

A more detailed comparison with iNakar. Piran & Granotl 
(2002) was performed and the results are listed in Table [2 
which includes the ratio of the numbers of observable or- 
phan afterglows to total observable afterglows. We choose the 
same SFR model (see Eq. (13) in Nakar, Piran & Granot 
2002), and the same model A and model B (i.e., the an- 
gle in Lorentz transformation 9 — 9 b s for model A, and 
9 = max(0, 6 b s — 6j) for model B), and neglect sideways ex- 
pansion. We conclude that the ratios for model B are all higher 
than thos e for model A. Our ra t ios are somewhat lower than 
those in Nak ar. Piran & Granotl (|2002J). This is because, for 
a fixed half-opening angle, the on-axis afterglow is brighter 
than in the sideways expansion case and the flux density of 
orphan afterglows does not increase very much, being due to 
the Doppler effect (see Eq. 0). 



The form of the SFR is 



10 0.75 Z 2 < 2pcak , 

b < 10°- 75z p- z > 2pcak , 



Opgi 




Fig. 3. The estimated number of orphan afterglows in a snap- 
shot for the whole sky, as a function of the limiting flux den- 
sity of detectors. The solid line represents our standard param- 
eterized result, with Ej = 1 x 10 51 erg, n = 1cm -3 , p = 
2.2, a 2 = 1.8, e e = 0.1, e B = 0.01, v = 4.55 x 10 14 Hz, and 
a power-law distribution of half-opening angles of the jets with 
^,min = 0.01, 0j, roax = 1. The thick dashed line and the thin 
dashed line are respectively the can onical line and optimistic 
line in lNakar. Piran & Granotl J2002I) . who assumed a laterally 
spreading jet. The jets have fixed initial half-opening angles 
9j =0. 1 and 0.05. The dotted line denotes the same parameters 
as the canonical line in lNakar. Piran & Granotl J2002I) . but the 
assumption of no sideways expansion is used. The dot-dashed 
line is the same as the d otted one except for the SFR model in 
|Porciani&Mad^(l200l). 



Recently, iRvkoff et alJ J2005I) performed a search for or- 
phan afterglows, but none has been detected. They gave an 
upper limit for the observed rate r? nia y < l.Qdeg^yr" 1 
by using the method suggested by Becker et al. d2004 : 



where B is the normalization factor. 



i] = N/((e)E) events deg~ yr~\ where N is the number 
of detected orphans, E the exposure, an d (e) the efficiency . 
Assuming a 30-minute exposure time as in lRvkoff et al.l J2005). 
we obtain the exposure E ~ 2.35 deg 2 yr for one whole 
sky survey. If the theoretical efficiency (e) is assumed to be 
unity, the observed rate is then 7/ ~ 1.3 x 10 _2 deg _2 yr _1 
and the detectability N is extrapolated to the 20th magnitude 
(N ~ 0.03) in the standard model. This is well below the upper 
limit 1.9deg _2 yr _1 estimated by Rvkof fetalJ d2005l) for the 
survey with limiting magnitude 20. 

We also calculated the detectability for different parame- 
ters to show their effects. Fig. 0] shows the detectability with 
different total kinetic energy Ej . The solid line is the standard 
one, the same as the solid line in Fig. [5] If the value of Ej in- 
creases by one order of magnitude, the detectability increases 
by a factor of about 12. It is reasonable that N is greater for a 
higher kinetic energy. Generally speaking, the detectability is 
sensitive to the Ej . 

Fig-0shows the detectability for different p and e e , with the 
solid line still the standard one. Remarkably, these parameters 
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M lim (Magnilude) 



7. 




F vj . m (ergcm-VH Z -') 



Fig. 4. Same as Fig. [5] with the standard parameters but for dif- 
ferent total kinetic energies. From top to bottom, the EjS are 
1.0 x 10 52 erg, 5.0 x 10 51 erg, 1.0 x 10 51 erg, 5.0 x 10 50 erg, and 
1.0 x 10 50 erg, while the solid line is the standard one. 



M ljm (Magnilude) 




F >Jim (ergcmVH z -') 

Fig. 5. Same as Fig. [3] The parameters are the same as the stan- 
dard ones (solid line) but for p = 2.4 (the dot-dash-dashed 
line), p = 2.1 (the dot-dot-dashed line), p = 2.05 (the dashed 
line), and e e — 0.3 (the dotted line). 



have minor effects on the results. Note that the p and a-i are 
independent paramet ers here, unlike the si deways expansion 
case, where c*2 = p JGranot & Kumarll2003h . so the variation 
in p does not change the temporal index a^. The parameter p 
influences the flux density at time tj, which can be seen in Eqs. 
and (|4}. The flux density is approximately proportional to 
e e , and the detectability is somewhat sensitive to e e . From those 
two equations, we can find that the other parameters n and €b 
do not significantly influence the detectability of optical orphan 
afterglows. 



4. Discussions 

We have calculated the number of orphan optical ( R-band) af- 
terglows that can be detected in an ideal survey of the whole sky 
with different limiting magnitudes. We considered jets without 
sideways expansion during the afterglow phase. This leads to 
a flatter light curve after the jet-break time tj, compared with 
the sideways-expansion case. Thus, orphan afterglows can per- 
sist for a longer time (above a certain flux density), and more 
expected orphans can be detected. The distribution of half- 
opening angles of jets was also considered, which suppresses 
the number of optical orphans compared with the case that 
all jets have one single half-opening angle 0j = 0.1. When 
combining these two effect s, the detectabi l ity is l ess than the 
canonical results of Na kar. Piran & Granot (2002), who con- 
sidered the case in which all jets have an initial half-opening 
angle 6j =0.1 with sideways expansion. 

From Figs. 0-(E}, we can conclude that the main fac- 
tors for the detectability are the total kinetic energy Ej, half- 
opening angle of jet 9j and te mporal index a? As Ej is consid- 
ered to be a standard energy ( Fr ail et al.ll200ll) . our results can 
be used to constrain the value of the Ej by detecting orphans. 
If Ej and a2 are determined accurately by other methods, the 
distribution function of the jet's half-opening angles can be 
determined well by observations of orphan afterglows. 

However, our estimation is simplified in several aspects. 
First, the parameters are undetermined and diverse in differ- 
ent bursts. For a variable parameter, we should know its dis- 
tribution. But this is very difficult. For example, the circum- 
burst environment seems to be an ISM, a wind, or another 
density-profile media, but these media cannot be determined 
clearly in well-observed GRBs. For simplicity, we choose 
the ISM with number density n = 1cm -3 . Second, or- 
phan afterglows may have propertie s similar to other phe- 
nomena: e.g., failed gamma-ray bu rst dHuang. Dai & Luf2 002: 
lRhoads ]|2003t iH uang et al.1l2005l) and "on-axis orphan after- 
glow" (Nakar, & Piran 2003). Third, we use the model with a 
power-law distribution of half-opening angles of uniform jets. 
However, ther e are some other st ructured jet models that cannot 
be ruled out dOai & GoijEoOll iRossi. Eazzati fc Reesll2002t 
Zhang & Mezsaros 2002). These models also affect the de- 
tectability. Fourth, we assume that the GRB rate is proportional 
to the SFR, so our results are SFR-dependent. Fifth, when ob- 
serving, one should distinguish orphan aftergl ows from other 
transients carefully JLevinson et alJl2002t iBecker et alJ feoO^ 
iGal-Yam et aklEoOfil) . Finally, the dust grains within the jet's 
opening solid angle may be evaporated by the prompt UV/X- 
ray photons, and the dust is possibly opaque in optical/UV 
bands outside the jet cone, so an optical orphan afterglow may 
be generally suppressed. 
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